Possible cooling by resonant Fowler-Nordheim emission 



Alexander N. Korotkov and Konstantin K. Likharev 
Department of Physics, State University of New York at Stony Brook, Stony Brook, NY 1 1794-3800 

(February 1, 2008) 



A new method of electronic refrigeration based on reso- 
nant Fowler-Nordheim emission is proposed and analyzed. In 
this method, a bulk emitter is covered with a-few-nm-thick 
film of a widegap semiconductor, creating an intermediate 
step between electron energies in the emitter and in vac- 
uum. An external electric field tilts this potential profile, 
forming a quantum well, and hence 2D electron subbands at 
the semiconductor-vacuum boundary. Alignment of the low- 
est subband with the energy levels of the hottest electrons of 
the emitter (a few ksT above its Fermi level) leads to a reso- 
nant, selective emission of these electrons, providing emitter 
cooling. Calculations show that cooling power as high as 10 4 
W/cm 2 (at 300 K), and temperatures down to 10 K may be 
achieved using this effect. 



The idea of using thermionic transport of electrons 
over an energy barrier for cooling has been repeatedly 
discussed in the literature (see, e.g., If the bar- 

rier height is a few times the thermal spread fcgT, the 
thermionic current may be quite substantial, with only 
the hot fraction of electrons being removed from the con- 
ductor. Unfortunately, the practical implementation of 
this idea runs into problems. 

A barrier of the necessary height (^200 meV for 300 K, 
and proportionally lower for lower temperatures) may be 
readily implemented in solid state structures, in partic- 
ular using composite semiconductors with the necessary 
conduction band edge offset. However, even if the bar- 
riers are relatively thick, the back flow of heat to the 
cooled conductor is prohibitively high |p]- ^|; multilayer 
structures proposed to overcome this effect [|] seem very 
complex and promise only a little cooling power. (Only 
at very low, millikelvin, temperatures where electron- 
phonon coupling is very weak, has efficient cooling been 
demonstrated using thermionic transfer over the super- 
conductor energy gap [@||.) 

Even a very narrow (submicron) vacuum gap can effec- 
tively quench the back heat flow, reducing it to radiation- 
limited levels of the order of 0.1 W/cm 2 (at 300 K). 
Unfortunately, in this case the energy barrier height is 
determined by the conductor workfunction which is too 
high (^> ksT) for most materials. A natural way to en- 
force electron transfer through a relatively high barrier 
is to apply a strong electric field (~ 10 MV/cm), in- 
ducing Fowler-Nordheim tunneling through the barrier. 
However, in typical situations the tunneling through the 
initially uniform barrier pulls out electrons within a rel- 
atively broad energy range that results in heating rather 



than in cooling (the "Nottingham effect" |IJ). 

We propose to limit the energy range of transferred 
electrons using resonant tunneling in a simple structure 
(Fig. |l|) where the bulk emitter (a metal or a heavily 
doped semiconductor) is covered with a thin (a-few-nm) 
layer of a widegap semiconductor. While at zero volt- 
age the electron potential energy profile of this structure 
has two steps (Fig. |l](a)), its tilting by the applied electric 
field creates a triangular-shape potential well at the semi- 
conductor film surface (Fig. 0(b)). Quantization of the 
energy of electron motion perpendicular to the surface 
results in discrete levels (subbands for the full energy) 
localized at the surface. If the electric field aligns the 
lowest subband with energy levels of the hottest electrons 
in the emitter (a few ksT above the Fermi level), reso- 
nant tunneling of these electrons to vacuum may lead to 
very efficient heat removal, and hence to emitter cooling 

Our proposal hinges on several ideas put forward ear- 
lier. First of all, numerous experiments indicate that 
Fowler-Nordheim emission from bulk cathodes is fre- 
quently enhanced by resonant tunneling through local- 
ized surface states arising from its unintentional contam- 
ination (see, e.g., Refs. ]8|||). Cooling of the nanoclusters 
using this effect was proposed in Ref. [[l0| . (Cooling of 2D 
electron gas based on the resonant tunneling through spe- 
cially fabricated quantum dots was proposed even earlier 
]ll[). However, to extend cooling to macroscopic objects, 
a large number of surface nanoparticles should be used 
in a single device. In this case, unavoidable spread of 
the size and shape of these particles would result in fluc- 
tuations of the resonant level positions, preventing their 
proper alignment with the Fermi level of the emitter, un- 
less nanoscale fabrication with atomic precision is used - 
a very distant prospect indeed. In contrast, our sugges- 
tion involves only planar structures and does not require 
nanofabrication. 

Concerning planar structures, Fowler-Nordheim tun- 
neling via the resonant subbands was predicted long ago 
|l2| and then observed in several systems Q-Q. The 
emission via resonant subbands at the outer surface of 
a widegap semiconductor in a strong electric field was 
predicted in Ref. pgfl . The Fowler-Nordheim emission 
coupled to the electron resonance in the vacuum gap was 
considered in Ref. [^0| (the implementation of this effect 
would, however, require an impractically fine uniformity 
of the gap). However, the possibility of heat removal was 
not mentioned in any of these publications. 

The objective of this work was a quantitative analy- 
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sis of the cooling effect in the system shown in Fig. [j], 
using a simple but natural model. First of all, we as- 
sume the interfaces to be perfectly plane. In this case 
the electron motion in the direction of tunneling (x-axis) 
and in the perpendicular direction (along the interfaces) 
are separated. Neglecting band bending due to quantum 
well charging (i.e. assuming its shape to be triangular) 
we have the following well known result for the electron 
eigenenergies (see, e.g., Bp): 



£ = £ x + £x = U — eEd + £ n , 
£ n = {-a n )(e 2 E 2 h 2 /2m) 1/3 



(1) 
(2) 



where all energies are relative to the emitter Fermi level, 
U is the initial energy step (Fig. |]), E is the electric field 
in the film, d is the film thickness, £± — h 2 k 2 _/2m, m 
is the electron effective mass in the conduction band of 
the film, and a n is the sequence of Airy function zeros: 
a = -2.34, ax = -4.09, . . . , a„ -> -[37r(n + 3/4)/2] 2 / 3 . 

In absence of energy relaxation, the level filling prob- 
ability p — p n (£±) may be found from the station- 
ary solution of the usual master equation, giving p = 
Jlhlilh +7k)j where / = /(£) is the Fermi distribution 
of the emitter electrons, and 7l and 7^ are the rates 
of electron escape from the quantum well into conduc- 
tor and into vacuum, respectively. These rates may be 
calculated as "/l,r = vDl^r, where v is the "attempt fre- 
quency", v = [2 / dxjvix)]- 1 = £ n /2h\a n \ 3 / 2 , and D LM 
are transparencies of the left and right triangular barriers 
[Fig. |l|(b)]. Within the WKB approximation (neglecting 
the image charge effects), 



\nD L = -^{eEd-£ n f'\ 
SeEn 



lnD R = -^^(<S>-U-£ n -A£f/ 2 . 
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Here the shift A£ = {Ji 2 k\j2)(mr 



(3) 
(4) 



is due to the 



possible difference between m and the electron mass mo 
in vacuum, <I> is the work function of the bulk emitter, 
and Eq is the electric field in vacuum. The relation be- 
tween this field and E includes the 2D charge density a 
of the electrons accumulated in the well: 



e E = ee Q E + a, 



(5) 



(e is the dielectric constant of the semiconductor film). 
The charge density <r, as well as the resonant current 
density j and thermal flow q, may be calculated as 



E 



ep, J 



E 

n.S 1 



e lap, q 



RP- 



(6) 



When the quantized level is above the emitter Fermi 
level, the typical spread of £± for the electrons in the 
subband is of the order of temperature T (here and below 
measured in energy units). Hence, assuming that the 



barriers are much higher than T, we can neglect the term 
A£ in Eq. (|J). Then integrating over £±, we get: 

j = epJ2 T\n{l + e- £ *' T ) lLlR /{ lL + lR ) 1 (7) 
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where Li2(z) = J^feLi z k /k 2 is the dilogarithm function 
and p = m/Trh 2 is the 2D density of states per unit area. 

Equations (0), (||) do not include the components of 
current (j') and heat flow (g') which are due to nonreso- 
nant, direct tunneling through the complete energy bar- 
rier. For this process, the barrier transparency may be 
calculated as D = DlD r . A standard WKB calculation 
yields: 
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, t = nT/£ , (10) 
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where l/£ = d{-hxD)/d£ =1/£q + l/£*, 

£^ = efiE/2{2m) 1 l 2 [U 1/2 - max(0, U - eEd) 1 ' 2 ], 
£§ = ehE /2(2m ) 1/2 (<i> - eEd) 1 12 , 
4(2m)V 2 



\a.Dk = — 



ZeEh 



[{U - £ x f /2 - max(0, U - £ x - eEd) 3/2 }, 



lnD« = - ^f^ max(0, $ - eEd - £ x f' 2 . 

(Notice that these formulas may be used even if £ x > 
U — eEd. The exact calculation in this case would give 
an extra factor on the order of unity, however, it can be 
neglected within the accuracy of WKB approximation.) 
The well-known factor p2|| tj sini shows that our approx- 
imation, based on the linear expansion of In D near the 
Fermi surface, can be used only at T < £ . (At higher 
temperatures the transport at large energies prevails.) 
At "low" temperatures we are discussing, the nonreso- 
nant tunneling always provides heating of the emitter, 
although it changes to cooling at the inversion tempera- 
ture T mv = £ /2. 

Figure || shows one of the results of our calculations 
using Eqs. The cooling power first increases ex- 

ponentially with the field, because the lowest subband is 
aligned with more and more populated hot electron lev- 
els, and then drops sharply as soon as the subband ap- 
proaches the Fermi level (at larger fields q becomes neg- 
ative). Just before this drop, the cooling power reaches 
a maximum, in this case, as high as 3 x 10 3 W/cm 2 at 
T = 300 K. 

The maximum values of g, as well as the correspond- 
ing values of j and q' , for several other parameter sets 
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are listed in Table I. From Fig. [2] one can see that the 
suitable range for the electric field E shrinks rapidly as 
emitter temperature goes down. Nevertheless, our model 
indicates that q + q' may be positive (i.e., cooling is still 
possible) for temperatures as low as 10 K. 

Let us discuss how realistic our model is. First of all, 
Eqs. (0) are valid only if the energy relaxation in the 
well is much slower than jl, 77?. We have also neglected 
the resonant subband broadening due to tunneling, but 
it was monitored through our calculations, so that the 
results are presented for only such parameters that the 
broadening is negligible in comparison with the impor- 
tant energy scale, T. One more possible source of devi- 
ations from our model is electron scattering in the well 
and during tunneling. However, these processes cannot 
lower the lowest subband deeper into the well, and so can 
hardly affect the process of hot electron extraction. Next, 
we have implicitly assumed that the Fermi energy of the 
bulk emitter is much larger than all the considered ener- 
gies. If this condition is not satisfied, the results would 
change, but not significantly. 

Despite the used assumptions, we expect that for 
smooth films the overall accuracy of our results for j, 
q, f, and q' for a given applied field is limited mainly by 
that of the WKB approximation |2^] and can be char- 
acterized by a numerical factor of the order of two or 
three. On the logarithmic scale at which we are working 
(Fig. ||) this is good accuracy indeed. Since the results 
show that the resonant emission cooling at temperatures 
above ~100 K may prevail over the nonresonant heating 
in a relatively broad range of electric field, and their ratio 
may be very high, we are confident that the net cooling 
of the emitter may be achieved. However, the estimate of 
the lowest achievable temperature (10 K) may be more 
vulnerable. 

The largest problem we see with the experimental im- 
plementation of resonant emission cooling is the neces- 
sary film uniformity. In fact, Table I shows that at 300 
K the effect is stable with respect to substantial (~ 20%) 
variations of d. However, to achieve cooling to 100 K, film 
thickness variations should not exceed ~ 4%, i.e., about 
one monolayer. (Thickness fluctuations require the elec- 
tric field to be decreased below the optimal value in order 
to be sure that we have not stepped into the heating re- 
gion on any considerable fraction of the emitter area.) 

While the proposed device potentially offers very large 
cooling power, its efficiency (more exactly, the "coeffi- 
cient of performance", COP) may be relatively low. For 
our system, COP may be presented as (q + q')/V(j +f), 
where V is the applied voltage. Even if the vacuum 
gap do is as small as 10 nm, the necessary voltage 
V = Ed + E d exceeds 10 V, giving COP about 10~ 3 at 
T = 300 K. To increase the COP, the electric field may be 
provided by a micromachined "grid" electrode very close 
to the surface, followed by another, much more distant 
grid at approximately the same voltage, and a collector 



at lower voltage, so that electrons are decelerated before 
the landing (see, e.g., p4|). 

To summarize, we have proposed a new method of elec- 
tronic refrigeration using the resonant Fowler-Nordheim 
tunneling in a fairly simple structure. If the experiment 
confirms our theory, this device may be invaluable for 
heat removal from electronic chips, as well as for the inte- 
gration of advanced low-temperature devices with room- 
temperature circuits. 

Useful discussions with D. V. Averin, H. Busta, and 
R. Tsu are gratefully acknowledged. 
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FIG. 1. The energy diagram of the proposed device: (a) 
in absence of bias and (b) at finite electric field. Resonant 
tunneling via quantized levels above the Fermi energy removes 
the hot fraction of electrons, thus cooling the emitter. 

FIG. 2. Solid lines: the resonant current density j, the 
corresponding cooling power density q, nonresonant current 
j' , and the corresponding heating power —q' as functions of 
the applied electric field E for $ —4 eV, U —1 eV, m — 0.5mo, 
m c = mo, e = 5, and d =2.5 nm, at T = 300 K. The dashed 
lines show the cooling power q at T = 100 K and T = 30 
K (the curve for —q' does not change significantly when the 
temperature is lowered). 
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TABLE I. Maximum cooling flow density q and heating 
density — q' , as well as the corresponding electric field E and 
resonant electric current density j, for several paratemer sets. 
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